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Abstract. The effect of the final state Coulomb interaction on particles produced in Pb+Pb collisions at 
158 A GeV/c has been investigated in the WA98 experiment through the study of the n~ /tt + and K~ /K + 
ratios measured as a function of mr - mo- While the ratio for kaons shows no significant mr dependence, 
the 7T~ /n + ratio is enhanced at small mr - mo values with an enhancement that increases with centrality. 
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A silicon pad detector located near the target is used to estimate the contribution of hyperon decays to 
the n~ /ir + ratio. The comparison of results with predictions of the RQMD model in which the Coulomb 
interaction has been incorporated allows to place constraints on the time of the pion freeze-out. 
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1 Introduction 

The distributions of negatively and positively charged pi- 
ons produced in heavy ion collisions exhibit significant dif- 
ferences at low transverse kinetic energy (or transverse 
mass), rriT — = \/v\ + m \ ~ m Tr- These differences 
are evident in the behaviour of the ratio of their yields 
R. K =ir~ /ir + . The study of charged particle production in 
Au+Au collisions at SIS (1 A GeV) [T] and AGS (10.8 A 
GeV) |2I3I4 | shows that for central collisions rises at 
low rriT at both collision energies. At the same time, for 
peripheral collisions at the AGS, does not depend on 
ttit- At 1 A GeV R^ is about 2.9 for tot — to w near zero 
and its integrated value is i? ff =1.9±0.1. The large inte- 
grated value at 1 A GeV can be explained by isobar decays 
and reflects the N/Z asymmetry of the colliding system. 
The very sharp rise at low tot on the other hand points 
toward significant Coulomb interactions between charged 
pions and the remaining nuclear charge distribution (ac- 
celeration of 7r + and deceleration of 7r _ ). At 10.8 A GeV 
and higher energy the isobar contribution to R n becomes 
insignificant and it is expected that mainly the Coulomb 
interaction distorts the R^ ratio. The NA44 experiment 
has published data on R^ measured in central Pb+Pb col- 
lisions at 158 A GeV A prominent enhancement of R^ 
at small tot was observed. NA44 has also measured R^ in 
S+S and S+Pb collisions at 200 A GeV/c, where the ratio 
remains constant at m-r — < 0.5 GeV/c 2 |2J. At high 
energies (AGS and SPS) the contribution from strange 
particle decay can significantly affect the measured pion 
distribution. The hyperon contribution is influenced by 
the detector acceptance and the pion vertex reconstruc- 
tion properties. The excess of hyperon over anti-hyperon 
yield in heavy ion collisions also leads to an enhancement 
of R-x at low transverse kinetic energy. 

Several dedicated models 6.7 8.9 l fTllllj describe the 
pion ratio detected by the NA44 collaboration in Pb+Pb 
collisions through the electromagnetic interaction induced 
by the large amount of charge of the participating pro- 
tons. The authors of the dynamical model of Ref. |S] ar- 
gue that the R v enhancement at small momenta depends 
mainly on the time of freeze-out and extract a freeze- 
out time of about 7 fm/c. They also predict the corre- 
sponding ratio for kaons, Rk = K~ / K + , to be smaller 
by a factor of to^ /"V . However, this overestimates the 
NA44 kaon measurement. This discrepancy is explained 
as due to the much larger absorption of K~ on nucle- 
ons {K~ + N — > A + 7r) than for K + . For a detector 
located at mid-rapidity this model predicts a small reduc- 
tion of the R v enhancement with respect to the NA44 
result measured near rapidity ycMS = 1 in the center of 
mass system. On the other hand, in another model |llj the 
enhancement of the ratio strongly increases towards mid- 
rapidity, reaching the value of 2.5. In ref. [5| it is concluded 
that the influence of the Coulomb force should be com- 



puted within an event simulator which takes into account 
a more accurate description of the space-time evolution of 
the collision. 

Initial results from the SPS experiment WA98 on the 
measurement of the ir~ /ir + ratio in central Pb+Pb colli- 
sions were reported in Ref. ^2- The enhancement of R v 
at rriT — <50 MeV/c 2 measured with the two inde- 
pendent tracking arms of WA98 located nearer to mid- 
rapidity than the NA44 spectrometer was reported to be 
significantly smaller than that measured by NA44 5 . 

In this paper we present the tt~ /tt + ratio as a func- 
tion of centrality for 158 A GeV/c Pb+Pb collisions. The 
background pion contribution from hyperon decays is de- 
termined and removed by use of vertex information from 
a silicon detector located near to the target to obtain the 
ratio of pions that originate in the target. The yields of 
electrons and positrons are analysed at low momentum 
where they can be clearly identified through time-of-flight. 
Their yields provide a check of the normalization of the 
opposite sign yields at the lowest transverse momenta. 
They also confirm the validity of estimates of the elec- 
tron and positron contaminations to the measured pion 
distributions at higher momenta. The K~ / K + ra.tio is also 
presented for central collisions. The results are compared 
with predictions of the RQMD 2.4 model ^Sj in which fi- 
nal state Coulomb interactions have been added. Within 
the uncertainties of the hyperon decay contribution the 
results are found to be compatible with predictions of the 
RQMD model in which the average pion freeze-out time 
is 15 fm/c for central collisions. Possible modifications of 
the final state predicted by RQMD, necessary to improve 
agreement with the measured tt~ /tt + ratios, are investi- 
gated and discussed. 



2 Experimental Setup and Data Analysis 

The CERN SPS experiment WA98 Q3| combined photon 
and hadron spectrometers with other large acceptance de- 
tectors that measured a number of global variables on an 
event-by-event basis. Fig. ^ shows the experimental lay- 
out during the 158 A GeV 208 Pb beam run period in 1996 
with a 0.239 g/cm 2 Pb target. 

The Mid- Rapidity Calorimeter, MIRAC [III, measured 
the collision transverse energy (Et) in the pseudorapid- 
ity interval 3.5 < rj < 5.5. The minimum bias trigger 
required Et 5 GeV and a valid signal from the beam 
counters. The measured da/dEx distribution is used for 
the calculation of the collision centrality. The minimum 
bias cross section for the run period used in this analysis 
is a m b=6451 mb with an overall systematic error of less 
than 10%. 

Two thin circular silicon wafer detectors were posi- 
tioned at small distance downstream from the target, where 
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Fig. 1. The WA98 experimental setup. 



charged particle trajectories were not affected by the mag- 
netic field. A Silicon Pad Multiplicity Detector, SPMD [Tfj| . 
was located at a downstream distance of 32.85 cm. It was 
segmented into 184 azimuthal sectors and 22 pseudorapid- 
ity rings in the range 2.35 < r\ < 3.75, and maintained pad 
occupancy below 20%. The SPMD measured the energy 
deposited in each pad. A Silicon Drift Detector, SDD 17 , 
was located at a distances of 12.5 cm from the target. 
Its position resolution was 25 and 35 /im in the azimuthal 
and radial coordinates, respectively. The acceptance of the 
two detectors overlapped in the pseudorapidity interval 
2.35 < rj < 3.4. The two detectors were used for precise re- 
construction of the vertex position of the Pb+Pb primary 
collision. For that purpose straight lines going through 
each pair of hits in the SDD and SPMD were projected 
onto the target plane. The one dimensional spatial reso- 
lution was found to be better than 0.3 mm perpendicular 
to the beam line. 

The WA98 experiment comprised two charged parti- 
cle spectrometer arms located on the right (first arm) and 
the left (second arm, not used in the present analysis) fac- 
ing downstream from the target, beyond a dipole magnet 
(Goliath) with 1.6 Tm bending power in the horizontal 
plane. The first tracking arm consisted of six Multistep 



Avalanche Chambers, MSAC, read out by CCD cameras 
equipped with two image intcnsifiers |18| . Each pixel of a 
CCD viewed a 3.1x3.1 mm 2 area of a chamber. In addi- 
tion, a 4x1.9 m 2 Time-of-Flight wall positioned behind 
the chambers at a distance of 16.5 m from the target 
allowed for particle identification with a time resolution 
better than 120 ps. The TOF detector consisted of 480 
scintillator counters with area 3.3x48.5 cm 2 and thick- 
ness 2 cm arranged in four rows |19|. Each counter was 
equipped with two PMTs, one at each end. The position 
of a hit along a scintillator bar was evaluated using two 
methods: the first used the time difference of signals ob- 
tained from the two ends of the scintillator bar, and the 
second used the ratio of their amplitudes. Both methods 
provide a spatial resolution of the order of 2.5 cm. 

Tracks were selected which traverse all six MSAC cham- 
bers, have detected hits in at least four of them, and were 
associated with a TOF hit within a 6x6 cm 2 window. An 
additional momentum dependent time of flight cut was 
used to exclude misidentified pion, kaon, and proton tracks 
at a level better than 1%. The (tot — ni r ) -rapidity ac- 
ceptance for pions generated by Monte Carlo is shown in 
Fig. |21 The center of mass system rapidity is 2.9. The mo- 
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mentum resolution of the spectrometer was Apjp = 0.005 
at p = 1.5 GeV/c. 




Fig. 2. Acceptance of the first tracking arm in the (mT—m n , y) 
plane for pions. 

Positive and negative particles were detected in two 
sets of runs with opposite direction of the magnetic field. 
Data samples consisting of 0.58xl0 6 ir~ and 0.22xl0 6 7r+ 
tracks satisfying all quality requirements were used for the 
ratio measurement. 

Efficiency matrices were measured for a grid of 8x8 
cells across each MSAC chamber and for several run peri- 
ods. The overall tracking arm efficiency measured for neg- 
ative particles was 2% less than that for positive particles 
due to chamber instabilities, however, there was no signif- 
icant difference in the uit dependence of the efficiencies. 
The detailed analysis of the 7r~ and K~ measurement in 
the first tracking arm has been presented in a separate 
paper [2T3| . 

3 Method of Simulation 

Simulations have been used in this analysis to estimate the 
different contributions to the measured and Rk ratios, 
and to compare with the measured results. A detailed de- 
scription of the primary collisions, of the Coulomb final 
state interactions, and of the particle propagation through 
the experimental setup has been simulated. 

The primary Pb+Pb collisions were simulated with the 
RQMD 2.4 model rT3j event generator which reproduces 
well the light strange baryon yield measured by the WA97 
experiment |21| . but does not include Coulomb interac- 
tions. In RQMD the particles are treated semiclassically 
over the complete evolution of the system. This allows the 
calculation of the Coulomb forces between each particle 
pair as an "afterburner" applied to the RQMD output. 
The electromagnetic interaction of a particle i with an- 
other charged particle j was calculated as described in 
Ref. j23 by 



where j runs over all charged particles of the system. Here 
particle i has mass m,, charge e, and 4- velocity Ui = 
( 1 , v) / vl — v 2 . The tensor FK" is given by 

where X x is the relative 4-distance between particles i and 
j . The forces are calculated at each time step of the system 
evolution and the positions and momenta of the particles 
are updated with the effect of the Coulomb forces only for 
those particles which have undergone their last collision. 
RQMD events with Coulomb interactions added were used 
as input to the detector simulation. Delta electrons pro- 
duced by Pb ions in the target material before the collision 
were also simulated through the detector response. 

The GEANT 3.21 program |H] was used for the de- 
tector simulation. The target, beam line, Goliath magnet, 
SDD and SPMD detectors, and tracking arm have been 
described in detail. For the track reconstruction the same 
procedure as used in the treatment of real data was em- 
ployed. For every hit in the tracking arm all information 
on the particle and all its predecessors was stored for fur- 
ther analysis so that the history of every track could be 
traced back to the moment of the primary collision. 

4 Experimental Results 

4.1 Centrality selection 

The R„ ratio was studied for six centrality intervals de- 
fined by the transverse energy Et measured with the MIRAC. 
The value of b may be estimated from the measured Et 
according to the equation 

oo 

nb 2 = J {da/dE T )dE T , (3) 

Ex 

where da/dEr is the experimentally measured distribu- 
tion. Systematic comparisons of da/ dE^ with predictions 
of VENUS 4.12 and extraction of the number of par- 
ticipating nucleons or number of binary collisions are given 
in Ref. [23 • Table 14. II lists the centrality intervals used in 
this analysis. 

4.2 Particle ratios 

The two tracking arms of WA98 provide two independent 
measurements of the pion ratio and have been shown to 
be consistent fT^j . This study focuses on an analysis of the 
data obtained by the first tracking arm, whose acceptance 
for pions is shifted by about 0.5 of a unit of rapidity closer 
to the target fragmentation region relative to the accep- 
tance of the second arm. This provides more favourable 
conditions for tz/K separation in the region of large trans- 
verse mass where the normalization of the ratios is per- 
formed. At low pt it allows better rejection of electrons 
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Centrality 

interval 

number 


(GeV) 


(%) 


Impact 

parameter 

(fm) 


Number 
of parti- 
cipants 


1 


34.5-130 


36.8-70.2 


8.6-12.0 


34-136 


2 


130-180 


28.0-36.8 


7.5-8.6 


136-179 


3 


180-240 


19.6-28.0 


6.3-7.5 


179-231 


4 


240-325 


10.1-19.6 


4.5-6.3 


231-306 


5 


325-400 


3.58-10.1 


2.7-4.5 


306-364 


6 


>400 


<3.58 


<2.7 


>364 



+ ^1.4 
1.2 



Table 1. Centrality intervals defined according to the amount 
of Et measured in MIRAC. The number of participants was 
calculated with VENUS 4.12. 



through timc-of-fhght and a significantly larger number of 
tracks pass through the active area of the SPMD than is 
the case for the second arm. Fig. [31 shows the ratio 
normalized at transverse kinetic energy between 0.3 and 
0.76 GeV/c 2 for the six centrality selections of Table |4~T1 
It should be noted that the Coulomb interaction effect on 
the ratio is strongly reduced and the contribution from 
hyperon decays is absent in the normalization region. 

The e~/e + ratio can be used as a check of the absolute 
normalization of the ratio. For 0.8 < p < 1.2 GeV/c 
electron and positron tracks are well identified through 
time-of- flight. As imposed by the tracking arm acceptance 
these tracks have transverse momenta below 50 MeV/c. 
Since both electrons and positrons originate mainly from 
photon conversions in the target material their yields should 
be identical, except for a small additional electron con- 
tribution from the production of delta electrons. These 
purely electromagnetic processes are well understood, so 
that the ratio e~/e + is expected to be reproduced in sim- 
ulation with an accuracy of order 1%. The e~/e + ratio is 
presented in Fig. 01 as a function of centrality and com- 
pared to predictions. For each centrality the e _ /e + ratio 
is normalized with the same normalization used for the 
corresponding ir~/ir + ratio in Fig. Any difference be- 
tween data and simulation is a measure of the systematic 
error in the normalization of the measured R^ distribu- 
tions of Fig. Summed over all centrality bins the dif- 
ference between data and simulation of the e~/e + ratio 
is -0.8±1.2%, consistent with the expected accuracy of 
the comparison. This result indicates that the systematic 
error on the absolute normalization of R^ introduced by 
normalization of the ratio in the region 0.3 < m-r — < 
0.76 GeV/c 2 is small, and that the normalization error is 
dominated by the statistical error of the data used for the 
normalization (shown by the dotted lines in Fig. 12). 

The ratio for kaons Rk measured under the same ex- 
perimental conditions is presented in Fig. 03 for the 10% 
most central collisions. The kaon ratio shows no signif- 
icant dependence on transverse kinetic energy. The lack 
of enhancement at low transverse mass is consistent with 
RQMD simulations with Coulomb final state interactions. 
The acceptance for the lowest transverse mass is located 
near kaon rapidity y=1.8 in the laboratory system. For 
m,T — nriK > 0.2 GeV/c 2 the average rapidity is about 
tf=2.15. 



+ ^1.4 
1.2 
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1.2 
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• Data 




0.2 0.4 0.6 0.2 0.4 0.6 

m T " m 7l 

( GeV/c ) 
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( GeV/c ) 



Fig. 3. The ratio R n vs tot — m v for several Pb+Pb collision 
centrality intervals defined by Et measured in the MIRAC. 
Data (filled circles) compared with predictions from RQMD 
including the Coulomb interaction (open circles). The ratios 
are arbitrarily normalized to 1.0 in the region 0.3 < niT-m, < 
0.76 GeV/c 2 . The statistical errors are shown by the error bars 
on the points. The dotted lines indicate the statistical error on 
the normalization of the data. The RQMD normalization errors 
are about half as large. 




400 
E T (GeV) 



Fig. 4. MIRAC Et dependence of the e~/e + ratio. The nor- 
malization is the same as for the tv~ /iy + ratio. Data are shown 
by filled circles. The hatched band indicates the statistical un- 
certainty (MeaniRMS) of the RQMD prediction. 
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Fig. 5. The kaon ratio for the 10% most central Pb+Pb colli- 
sions. Data (filled circles) are compared with predictions from 
the RQMD model with Coulomb interactions included (open 
circles). The ratios are arbitrarily normalized to 1.0 in the re- 
gion 0.3 < rriT — m,K < 0.6 GeV/c 2 . Errors are as described 
for Fig. OH 



4.3 Hyperon contribution 

Pions from hyperon decays have a much steeper trans- 
verse momentum spectrum than directly produced pions. 
Since strange hyperons are more abundantly produced 
than anti-strange hyperons, the tt~ /ir + ratio will be en- 
hanced at low rriT due to this difference. In order to study 
the Coulomb effect on the directly produced pions, it is 
necessary to estimate and remove the hyperon contribu- 
tions. 

The Silicon Pad Multiplicity Detector has been used 
to estimate the yield of identified tracks that do not track 
back to the target vertex, and thus may be attributed 
to products of strange or anti-strange hadron decays. All 
pads of the SPMD that overlap with a 4x4 mm 2 win- 
dow, centered at the point where the extrapolation of the 
track towards the target intersects the SPMD plane, were 
selected. The summed energy loss in the pad window is 
denoted as S, and given in units of the energy loss of 
a minimum ionizing particle in the silicon wafer of the 
SPMD. S therefore gives an estimate of the number of 
minimum ionizing particles traversing the SPMD within 
the pad window. Simulation indicates that the mean value 
over a large set of tracks for pion tracks originating in the 
target is (S) = 0.89, while for tracks from hyperon and 
Kg decays (5)=0.17 and 0.20, respectively. This is be- 
cause tracks that originate from hyperon or Kg decays ei- 
ther traverse the SPMD far from the pad indicated by the 
extrapolation, or do not traverse it at all (produced down- 
stream of the SPMD). Thus, tracks from hyperon decays 
will be suppressed by roughly a factor of 0.89/0.17 « 5 
with the SPMD hit requirement. 

In the case of high pad occupancy the energy loss 
Si measured in the SPMD window associated with the 
track contains a significant contribution 5*2 from spuri- 
ous tracks. The value of S2 is estimated by the measured 
energy loss in a window of the same size centered at a lo- 



cation far from the extrapolated track location (in partic- 
ular, at a location rotated by 12° in azimuth). The RQMD 
simulation indicates that (S") = (Si — Sa) provides a good 
estimate of (S) for each kind of track. 

For tracks whose extrapolation traverses the SPMD 
the quantities R g = ir~ /n + and Rspmd = S'~ / S l+ are 
defined, where S'~ and S' + are sums of S\ — S2 values for 
negative and positive pion tracks identified by the tracking 
arm. In the case of equal hyperon fractions of the pion 
yield for both charges one expects Rspmd ~ Rg- The 
behaviour of R g and Rspmd as a function of transverse 
kinetic energy for the 10% most central collisions is shown 
in Fig. El Both distributions have the same normalization 
chosen to give R g = 1 at mj — > 0.3 GeV/c 2 . 

Since the contribution to S from hyperon decay tracks 
is smaller than for tracks originating in the target, a devi- 
ation of Rspmd from R g is an indication of different hy- 
peron fractions in the two charge samples. The observed 
lower value of Rspmd m comparison to R g for small ttit 
indicates a larger fraction of tracks from hyperon decays 
in the measured n~ distribution relative to the 7r + dis- 
tribution. To quantify this difference the quantity (R g — 
Rspmd)/ Rg is extracted and compared with RQMD sim- 
ulation as a function of uit — in Fig. Hand as a function 
of collision centrality in Fig.[S]for mT — m^ < 140 MeV/c 2 , 
where the excess becomes significant. Within errors it is 
seen that the RQMD simulation provides a good descrip- 
tion of the contribution of strange baryon decays to the 
pion ratio. 
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Fig. 6. Rspmd (triangles) and R g (circles) for the 10% most 
central events. See text for definitions of Rspmd and R g . Both 



distributions are normalized 
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The yield distributions of light strange and anti-strange 
baryons have been reported by the WA97 collaboration 
[22] and also shown to be successfully reproduced by RQMD 
calculations. Based on this fact, and on the comparisons 
of Figs.[7|and |3 the RQMD simulation results have been 
used to correct the measured pion ratios R^ shown in 
Fig. |21 for the hyperon decay contributions to obtain the 
corrected ratios R'^ of pions originating in the target. For 
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Fig. 7. (Rspmd — R g )/R g for the 10% most central events. 
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Fig. 9. The corrected charged pion ratio R'^ = n~ /% + for 
pions originating from the target, after removal of the hyperon 
decay and electron contamination contributions. Data (filled 
circles) are compared with predictions from RQMD including 
the Coulomb interaction (open circles). The ratios are arbi- 
trarily normalized to 1.0 in the region 0.3 < mr — m n < 
E T (GeV) Q 75 GeV/c 2 . Errors are as described for Fig. El 



Fig. 8. Centrality dependence of (Rspmd — R g )/R g for 
rriT — < 140 MeV/c 2 . Data are shown by filled cir- 
cles. The hatched band indicates the statistical uncertainty 
(MeaniRMS) of the RQMD prediction. 



this purpose, for each of the ir~ and tt + yields, the fraction 
of the pion yield from strange hadron decays has been esti- 
mated from RQMD simulation and removed from the total 
yield. Fig.FJJshows the corrected pion ratios R'^ = ir~ /ir + . 

At low momenta the e - and e + tracks have been sep- 
arated from pions through time-of-flight and their abun- 
dances are found to be in agreement with simulation. At 
high momenta e~ and e + tracks cannot be separated from 
pions and their relative yields have been estimated from 
simulation. The ratio R'^ of pions originating from the tar- 
get has also been corrected for e~ and e + misidentified as 
pions at high momenta. According to simulation this e~ 
and e + contamination weakens the enhancement of the 
pion ratio by about 2% for the lowest rop. The corrected 
ratios R'^ shown in Fig. El can be used for comparison of 
the 7r~/7r + ratio with predictions from models including 
the decay of the 77 and other short lived resonances. 



In Fig. the dependence of the corrected pion ratio 
R'^ on the collision centrality is shown for two intervals 
of pion transverse kinetic energy: below 40 MeV/c 2 and 
between 40 and 140 MeV/c 2 . The error due to normaliza- 
tion is included in the error bars. Predictions from RQMD 
calculations with Coulomb interactions included are also 
shown. Both the measured results and the RQMD predic- 
tions show a smooth increase of the pion ratio with in- 
creasing centrality. As expected, the ratio tends to a value 
of unity for peripheral collisions. However, RQMD is seen 
to overpredict the measured pion ratio with a discrepancy 
that grows with decreasing transverse mass. 



5 Modified Coulomb Calculations 

Several parameters of the collision, such as the total par- 
ticipant charge at central rapidity, its transverse expan- 
sion velocity, and the freeze-out time, affect the it~ /n + 
ratio |H] . In order to investigate the sensitivity of the pion 
ratio to the collision parameters, the Coulomb calcula- 
tions were repeated under alternative scenarios following 
minor modifications of the RQMD output. For this study 
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ratio for pions originating from the target for bit — rn^ < 
40 MeV/c 2 and 40 < m T -m T < 140 MeV/c 2 . Data (filled cir- 
cles) are compared with predictions from RQMD calculations 
(open circles). 



RQMD events from central Pb+Pb collisions with impact 
parameter b <4.5 fm were used. 

In the first modification, the RQMD baryon rapidity 
density was adjusted. Since the value of dN/dy for pro- 
tons around mid-rapidity for central Pb+Pb collisions as 
measured by the NA44 J2H| and NA49 |22j experiments 
is lower than the RQMD prediction by about 30%, the 
RQMD Coulomb interaction from baryons should be re- 
duced accordingly. The corresponding reduction of the 
baryon Coulomb source was implemented in the calcu- 
lation by weighting the Coulomb charge of all baryons 
during all stages of the collision by the ratio of measured 
to predicted proton dN/dy distributions. This baryon ra- 
pidity density correction (BRDC) reduces the ratio at the 
lowest rriT — bin by 3% (Fig. 111(1 . Since this correction 
is based on a discrepancy between RQMD and experimen- 
tal measurements, further calculations use this modified 
baryon Coulomb charge rapidity density. 

To estimate the sensitivity of the pion ratio to the 
freeze-out time the application of the Coulomb interac- 
tion was simply delayed. The mean time of pion freeze- 
out at \ycMs\ < 1 in the central Pb+Pb RQMD events is 
15 fm/c. For demonstration, a 2 fm/c delay was added dur- 
ing which time the particles move freely after freeze-out 
on their final RQMD trajectories. After propagating freely 



1.3 



1.2 



1.1 



■+■ - RQMD and Coulomb 

• - baryon rapidity density correction (BRDC) 

* - BRDC and pion freeze -out delay 2 fm/c 



0.1 0.2 0.3 0.4 



0.5 



0.6 0.7 
( GeV/c 2 ) 



Fig. 11. The pion yield ratio -k~ at mid-rapidity (1 < 
\ycMs\)- Predictions from RQMD with final state Coulomb in- 
teraction unmodified (crosses), with BRDC (circles), and with 
BRDC and pion freeze-out delayed an additional 2 fm/c (tri- 
angles) . 



for 2 fm/c the pions are then subjected to the Coulomb 
interactions. At vtit — m w <40 MeV/c 2 the resulting ra- 
tio shown in Fig. 1111 by triangles is reduced with respect 
to the unmodified RQMD result by 5%. This total reduc- 
tion of the ratio is similar in magnitude to the difference 
between the data and the unmodified RQMD simulation 
results of Fig. Alternatively, increased collective flow 
in RQMD would result in a larger freeze-out volume and 
a reduced Coulomb interaction, with a similarly improved 
agreement with the measurements. 
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Fig. 12. Comparison of NA44 data (open circles) and WA98 
data (filled circles) for the 10% most central collisions. The ra- 
tio for pions originating in the target is also shown (triangles) . 



In Fig. 1121 the pion ratio measured in this experiment 
for the 10% most central Pb+Pb collisions is compared 
with that reported by the NA44 collaboration [Sj- Nei- 
ther RQMD calculations, nor the previously mentioned 
models [oTlll) . predict the observed difference in the two 
measurements. This could indicate an earlier pion freeze- 
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out in the center of mass rapidity region near y=l (NA44) 
relative to that near y=0 (WA98). 

The measured K~ / K + ratio (Fig. does not show 
any significant wit dependence. This is consistent with 
the data of the NA44 collaboration and with RQMD pre- 
dictions. The kaon ratio enhancement induced through the 
Coulomb interaction is expected to be less than that for 
pions by a factor equal to their mass ratio m^/niK [ZJ- 
This relation is satisfied with respect to the pion ratios 
presented here, but it is not for the NA44 data [5]. There- 
fore, while the Coulomb effect is consistent with the WA98 
kaon result, an additional mechanism, such as different ab- 
sorption of K + and K~ , is necessary to explain the NA44 
result. 



6 Discussion and Conclusions 

The 7r~ / 7r + ratio has been measured as a function of trans- 
verse mass and centrality in Pb+Pb collisions at 158 A 
GeV/ c. The hyperon decay contribution has been deduced 
and shown to result in an increase of the pion ratio at 
low rriT that is responsible for about half of the enhance- 
ment of tt~ relative to 7r + in central collisions, and an even 
larger fraction of the enhancement in peripheral collisions. 
The RQMD model was shown to agree with the measured 
hyperon decay contributions. The hyperon decay contri- 
butions have been removed from the measured ratios to 
extract the ratios for directly emitted pions. 

The low rriT enhancement of the tt~ /tt + ratio, after 
removal of the hyperon decay contribution, increases with 
collision centrality and tends to zero for the most pe- 
ripheral events (see Fig. EJ). This observation is consis- 
tent with the hypothesis that the enhancement is due 
to the Coulomb interaction induced by the net positive 
charge of the participant protons. The same behaviour is 
observed in the framework of RQMD model simulations 
that include final state Coulomb interactions. However, 
the RQMD simulation overpredicts the observed enhance- 
ment in the ratio by a factor of about 1.5. 

Investigations of the RQMD model predictions sug- 
gest that rather small modifications of the properties of 
the participant fireball are necessary to obtain good agree- 
ment with the data. For example, good agreement can be 
achieved after reduction of the baryon rapidity charge den- 
sity predicted by RQMD according to the NA44 and NA49 
proton measurements, together with a small increase of 
the mean pion freeze-out time (such as an increase from an 
average of 15 fm/c to 17 fm/c). According to RQMD pre- 
dictions, the transverse flow velocities of heavy particles 
are less than those of pions [2Hj- A slight increase of the 
baryon transverse flow velocity would have a similar effect 
as the delay of the freeze-out time. In both cases, the sys- 
tem disperses over a larger volume before freeze-out, which 
reduces the Coulomb field and, therefore, the ir~ /ir + ra- 
tio. Based on comparisons with the RQMD model cal- 
culations, the results suggest a relatively large value of 
the mean freeze-out time for pions (~15 fm/c) in central 
Pb+Pb collisions. 
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